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ABSTRACT 



We used the large binocular camera (LBC) mounted on the large binocular telescope (LBT) to observe the Lockman Hole in the U, B, 
and V bands. Our observations cover an area of 925 arcmin 2 . We reached depths of 26.7, 26.3, and 26.3 mag(AB) in the three bands, 
respectively, in terms of 50% source detection efficiency, making this survey the deepest U-band survey and one of the deepest B and 
V band surveys with respect to its covered area. We extracted a large number of sources (~ 89000), detected in all three bands and 
examined their surface density, comparing it with models of galaxy evolution. We find good agreement with previous claims of a steep 
faint-end slope of the luminosity functions, caused by late-type and irregular galaxies at z > 1.5. A population of dwarf star-forming 
galaxies at 1.5 < z < 2.5 is needed to explain the U-band number counts. We also find evidence of strong supernova feedback at high 
redshift. This survey is complementary to the r, i, and z Lockman Hole survey conducted with the Subaru telescope and provides the 
essential wavelength coverage to derive photometric redshifts and select different types of sources from the Lockman Hole for further 
study. 

Key words. Surveys - Galaxies: photometry 



1. Introduction 

The formation and evolution of cosmic structures, such as galax- 
ies, clusters, and the large-scale structure, are some of the most 
important issues in modern astrophysics. According to hierar- 
chical models, initial fluctuations of the dark matter mass den- 
sity develop to form galaxies, clusters, and the cosmic web. 
Such processes leave their footprints in different regimes of the 
electromagnetic spectrum, and assembling statistically signifi- 
cant samples of extragalactic objects at different wavelengths 
can give valuable information on the various processes involved 
in the evolution of the universe. 

A very valuable tool for constructing such samples is deep 
"blind" surveys, where a region in the sky with no bright sources 
is observed with a long integration time. Optical surveys are 



very important in this context, as they are able to provide the 
densest fields in terms of detected sources and serve as "an- 
chor points" for the multi-wavelength coverage. After a multi- 
wavelength coverage has been ac hieved, one could apply pho- 



tometric redsh i ft techniques (e.g. iBolzonel la. Miralles & P ello. 



2000} iBenftezt l2000t lllbert et al l I2009I) to examine the lumi 
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nosities of the various sources or select source samples for spec- 
troscopy 

Notable results have been reported in various fields of ex- 
tragalactic astrophysics using blind deep surveys. Combining 
imaging and spectroscopic surveys at different regimes of 
the spectrum, differe nt group s have been able to derive the 
star formation ( e.g. | Hopkinsl |2004|) and accretion histories 
(e.g. lUeda et all 2003b of the universe and examine their co- 
evolu tion dVollmer. Beckert & Daviesl l2008h ISomerville et all 
I2008I) . From optical imaging and photometry alone, one can use 
the information in the number count of the detected sources to 
test the geometry and evolutionary models of the universe. For 
example, Eucledian geometry would result in a constant slope of 
0.6 in the galaxy number counts with respect to their magnitudes, 
but t his has been ruled out from early re sults in this direction 
(e.g. iGardner. Cowie & Wainscoatl Il993l) . Measuring the num- 
ber counts in different wavebands, it is evident that simple ge- 
ometric models invoking a "deceleration parameter" (q) could 
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not give goo d fits and some kind of evolution has to be taken 
into account dMetcalfeetal.Lll991l) . This effect is more severe in 
blue colours in the form of excess counts at fainter magnitudes 
and it is widely known as the "faint blue gal axy problem " . With 
high resolution observations using the HST. lDriveret"ai] (119951) 
demonstrate that the sources responsible for the faint counts 
have late-type and irregular morpho logies; adding a popu lation 
of z — 2 dwarf star-forming galaxies dMetcalfe et all ll995) gives 
a reasonable fit to the blue number counts data. These galaxies 
contribute to the star formation at redshifts z > 1 and are merged 
or simply have evolved to non activity locally. 

Support for this scenario comes from the study of the (blue) 
lumin osity functio n s of d ifferent kinds of objects at different red- 
shifts. Illbert et al] d2005l) find that bluer luminosity functions 
show evidence of more rapid evolution with redshift than red- 
der ones, and later spectral ty pes and bluer colours seem to play 
a mo re important role in it (IZucca et all 120061; IWillmer et all 
2006). Howeve r, small evolution o f the disc population to z — 1 
is observed by Illbert et al] d2006[) , but the strong evolution of 
bulge-dominated sy stems could be attributed to a dwarf galaxy 
population (see also llm et alll200l"h . The study luminosity func- 
tion is limited to relatively bright objects as it is based on red- 
shifts. A number count distribution can probe fainter objects and 
give an approximation on the faint-end slope (Ba rro et al 1 12009b 
of the LF and help distinguish bet ween different resu lts (see 
comparisons in Illbert et alj|2005l and IZucca et al]|2006l) . In this 
paper we present deep U-B-V band observations of the Lockman 
Hole with the corresponding number counts to 27.5 mag(AB). 




3000 4000 5000 6000 7000 

MA) 

Fig. 1. Transmission curves of the filters used (U-spec - B-Bessel 
- V-Bessel) are shown by continuous lines. The dotted line rep- 
resents the U-Bessel filter available on the LBC. We have chosen 
the U-spec filter for the U-band on grounds of its better efficiency 
and more uniform spectral coverage. These curves represent the 
filters' responses without accounting for the detectors' responses 
or the atmosphere. The detectors' responses are slightly different 
for the two arms of the telescope which might have an effect in 
the V-band. 



2. The Lockman Hole multi-wavelength survey 

The Lockman Hole is a regi on with minimal galactic absorp- 
tion ( Nhi = 4.5 x 10 19 cm" 2 . iLockman. Jahoda & McCammonl 
1986) and the absolute minimum of infrared cirrus emission 
in the sky. Its position in the northern sky (a = 10 h 52 m 43 s , 
5 = 57°28'48") makes it an ideal location for deep surveys. 
Indeed it has a large multi-wavelength coverage spanning from 
X-rays to meter-wavelength radio. In X-rays it h a s bee n ob- 
served with the ROSAT satellite dHasinger et all Il998l) and 
more recently with XMM (Hasinger et al., 2001; Brunner et all 
2008), reaching a depth of 1.9 x 10" 16 ergcm- 2 s _1 in the 0.5- 
2.0keV band . In the ultra-violet it has been observed by GALEX 
(Mar tin et all 120051) as one of its deep fields, with the data be- 
ing publically available. In the near infra red (J and K bands) it is 
a part of the UKIDSS ultra deep survey dLawrence et all [20071) 
reaching K~23(AB). In infrared wavelengt hs it was observed by 
ISO using both ISOPHOT and ISOCAM ([Kawara et all l2004t 
Fadd a et all 120041; Rodighie ro et all 120041) and more recently 
there have been observation s with Spitzer-IRA C dHuang et all 
2004) and Spitzer MIPS dEgami et all 120081). The Lo ckman 
Hole is also part of the SWIRE survey dLonsdale et all 12003b . 
observed with both IRAC and MIPS and covering a much 
wider (but shallower) area. There have been a number of 
millimeter - sub-mm observati ons of the Lockma n Hole, 
namely with the JCMT-SCU BA dCoppin et all 120061). JCMT- 
AzTE C dScott et all l2006j). IRAM-MAMBO dGreve et all 
120041) , and CSO-Bolocam ([Laurent et all 120051) . In the radio 
regime, the Lockman Hole h as been observed w i th the VLA , 
both i n 5 and in 1.4 GHz dCiliegi et all 120031; llvison et all 
120021; iBiggs & IyisonL 120061) and with MERLIN in 1.4 GHz 
(Big gs & Ivisonll2008l). Finally, in met er-wavelengths it was tar- 
geted by the GMRT dGarn et all 120081) . 

In this work we present the results of an imaging campaign 
of the Lockman Hole in the optical. We have used the LBT to 



obtain deep U, B, and V images. The "red" part of the optical 
imaging campaign has been conducted with the Subaru telescope 
(in the r, i, and z bands) and will be presented by Szokoly et al. 
(in preparation). 



3. Observations 

The o bservations were made with the Large Binocular Camera 
(LBC, iGiallongo et all [2008) of the Large Binocular Telescope 
(LBT) on Mount Graham, Arizona. The LBT has two 8.4m 
mirrors on a common mount and both of them are equipped 
with a prime focus camera. Both LBCs contain four CCD chips 
with 2048x4608 pixels each. Three chips are aligned parallel to 
each other while the fourth is tilted by 90 degrees and located 
above them. This provides a 23x23 arcmin field of view with a 
sampling of 0.23 arcsec/pixel. The gaps between the CCDs are 
945 nm wide, which corresponds to 18 arcsec, thus a 5-point cir- 
cular dither pattern with a diameter of 30 arcsec was chosen to 
provide good coverage over the whole area. 

Both cameras have an 8 position filter wheel each, and to- 
gether a total of 13 filters are available, covering a range from 
the ultraviolet to the near-infrared. For the U-band Lockman 
Hole imaging we used the special LBT U-band filter (see Fig. 
[1} which has a more uniform coverage and better efficiency than 
the standard U-Bessel. For the other images we used the standard 
B-Bessel and V-Bessel filters. 

The Lockman Hole was observed in March, April and May 
2007 during science demonstration time (SDT; PI: E. Egami), 
when only the camera on the "blue" channel of the telescope 
was available, and in 2008 and 2009 during LBTB (German 
institutes') time (PI: G. Hasinger) with both cameras avail- 
able in "binocular" mode. We have chosen 2 different point- 
ings as centres of the image, corresponding to the VLA (a = 
10 h 52 m 08.8 s , 6 = 57°2 F 34". llvison et alj |2002) and the XMM 
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Table 1. Details of the observing runs (time is in minutes). The 
SDT runs were made using the blue channel only, whereas the 
LBTB runs used both telescopes. The effective exposure time is 
the time actually used for the final stacks. 



(a = 10 h 52 m 43 s , 6 = 57°28'48". lHasinger et alj|200ll) point- 
ings. These are separated by 8.6arcmin, so there is a large area 
of overlap, where our images have the highest sensitivity. During 
the science demonstration time the observing time was split in 
half between XMM and VLA exposures and during the LBTB 
time we concentrated on the XMM area. 

The total time spent on the Lockman Hole was 36.8 hours 
which are distributed among the various filters using both chan- 
nels (when available) as described in Tab. [3] The exposure time 
for each observation was 360 seconds initially, but it was reduced 
to 180 seconds for later observing runs (May 2008 onwards), af- 
ter discovering a large number of saturated sources and limited 
source tracking efficiency of the telescope for long exposures. 
The effective exposure time however is as observational prob- 
lems such as high altitude cirrus clouds or bad seeing diminish 
the quality of certain images which were not used for creating 
the final stacks. As seen from Tab. [3] the time efficiency of the 
three bands is in the order of 65%. 

We should note here that the V-band observations were taken 
using both the blue and the red arms of the telescope. The 
blue arm was used during SDT, and the red during LBTB time. 
Although the response curves of the two V-band filters are iden- 
tical, the quantum efficiencies of the detectors are slightly dif- 
ferent. For the analysis presented in this paper this effect is not 
significant and we merge the two (V/, and V r ) images to achieve 
greater depth. However in more detailed studies, one should treat 
the Vb and V r images separately. 



4. Data reduction 

For the reduction of the data we have primarily used tRAF rou- 
tines included in the mscred package, which is designed to re- 
duce mosaic data. 



4.1. Initial calibration 

Initial corrections to remove the "pedestal" level of each chip 
have been carried out using the overscan regions. We have 
found that the level of the corrections varies significantly (up 
to 5%) with the column of the chip and therefore we fitted an 8- 
order Legendre polynomial to it. Residual errors (possibly row- 
dependent) have been corrected using bias frames taken at the 
beginning and end of each night with zero integration time. 

Flat corrections have been made using sky flat frames taken 
at either dusk or dawn (or both) for each filter at each arm of 
the telescope separately. A master flat has been created for each 
night, filter, and arm. We divided the bias-corrected images with 
their respective flat fields and noticed that the outer edges of 
the fourth chip, corresponding to the outer edges of the field of 
view were extremely noisy, possibly due to poor illumination. 
We flagged them as bad regions. 

After flat-fielding, we corrected the images for bad pixels. 
These include columns of the CCD with non linear response or 
dust on the CCD surface. Bad pixel masks are created and the 
correction has been made by interpolating the values of neigh- 
boring pixels. Note that because of the absence of neighboring 
pixels on the edges of the field of view, these areas are not cor- 
rected and are simply not taken into account for the final stages 
when the images are stacked. 

Finally, there are some bright sources in the field which sat- 
urate the response of the CCD. In the most severe cases the flux 
is so high that the current affects the neighboring pixels, leav- 
ing "bleeding trails". This has a significant effect to the B and 
V images and therefore these regions have been identified and 
masked out. 



4.2. Astrometry 

Before dealing with the (arcsecond-scale) astrometric errors, we 
correct each image for an initial offset, in the order of several ar- 
seconds, caused by the telescope's pointing inacc uracy. We use 
the brightest (r < 19) sources from the USNO-A2 dMoned[l998l) 
catalogue to correct for this offset. This is done by simply updat- 
ing the wcs header of each file to match the coordinates of the 
catalogue stars. 

After having done that, we need to correct for the true as- 
trometric errors caused by the camera distortion. For this pur- 
pose we do not use the USNO catalogue of the brightest stars, 
as proper motions could have an effect in the solutions we de- 
rive. We use an astrometry corrected catalogue of the Lockman 
Hole, which includes sources brighter than V = 19. This is 
based on observations made with the Canada-France-Hawaii 
Telescope (CFHT , e.g. 



Wilson et al 112001b and the data reduc- 
99%~ 



tion details are in iKaiser etaL (1999). The absolute astrometry 
of this catalog ue is based on USNO-A2 which claimed accuracy 
is 0.25 arcsec dMoneA[l998h . 

To apply detailed astrometrical solutions we deal with each 
chip separately in order to avoid fitting for jumps between the 
chips. We first apply a distortion pattern which we empirically 
derived by correcting a random image and then fit a 4-order poly- 
nomial to each direction of each chip. The final rms scatter we 
get is in the order of 0.2 arcsec. An example of the astrometrical 
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Fig. 2. Astrometrical solutions applied to an exposure image of 
1 1 June 2007 (U-band). The image has been bias and flat cal- 
ibrated but no astrometrical solution has been applied. The ar- 
rows represent the shifts to the calibration sources exaggerated 8 
times. To compensate for the distortion we dealt with each chip 
separately and after applying an initial shift we fitted a 4th order 
polynomial to the corrections. The corrected images have an rms 
scatter of ~0.2 arcsec 

solutions applied (after correcting for the overall pointing offset) 
is given in Fig. [2] To measure the final astrometrical accuracy, we 
compare our LBT images with the USNO-A2 catalogue (see Fig . 
O and with others, such as USNO-B1, ATM dlrwin et al.Lll994l) . 
SWIRE-IRAC(3 .2um) Co nsdale et al.11 20031) and L-band VLA 
(Bigg s & Ivisonl |2006). Their positions typically agree within 
0.4 arcsec and the typical standard deviation is 0.45 arcsec, 
which is the value we assume to be our final astrometric ac- 
curacy. The relative astrometry of the U-B-V images based 
on the positions of bright (<24mag) and realtively compact 
(FWHM<1.5 arcsec) sources has a standard deviation of 0.066 ". 

4.3. Background subtraction 

The final step is to subtract the sky background. After having 
flattened the images and having corrected for field distortions 
(without preserving the flux of each pixel) the sky background is 
uniform within a good approximation. To subtract it we frag- 
ment the image to a grid constructed of 100x100 pixel wide 
meshes and smooth each mesh by a median filter with a 5x5 pixel 
kernel using sextractor (Bert in & ArnoutsH T996). We chose this 
method over fitting a function to the background because it gives 
better results in the vicinity of bright stars in the sense that it does 
not over-correct the background. 

After having subtracted the background, we re-project the 
four chips of each image to a common frame, applying the com- 
plex astrometrical solutions. By doing that we get rid of com- 
plex headers and multiple frame images. We use the same ref- 
erence image to re-project all images in all filters. Finally, after 



removing any bad images due to poor observing conditions or 
other problems, we stack all the images (weighted according to 
their exposure times) to produce the U, B, and V maps of the 
Lockman Hole. An example image (in the U band) and its cor- 
responding exposure map is shown in Fig. [4] The regions mark 
the deep VLA and XMM surveys with 10'radii, which are their 
typical widths. 

4.4. Flux calibration 

For the B and V filters we rely on flux-calibrated images of 
the Lockma n Hole taken with the Ca lar Alto Telescope (see 
iKaiser et all 119991; IWilson et all 120011) . We select point-like 
sources which are not saturated in any of the images and conduct 
aperture photometry. We compare the results and derive zero- 
point magnitudes for our final images. We do not find evidence 
for a gradient across the image. 

As there are no sources with known magnitudes in the U- 
spec band, we had to rely on U-Bessel standards to derive the 
zero-point offsets. We used the observations of June 11, 2007 
when standard stars are observed with both the U-Bessel and 
the U-spec filters. We have applied the same calibration (bias 
subtraction and fiat fielding) using the same bias and fiat-field 
images to all the target and standard star frames and did not per- 
form any astrometric corrections nor we combined the calibrated 
files. We measured the observed magnitudes of four standard 
stars without applying any zero-point offsets and found a dif- 
ference of ~ 0.67 mag. We attribute this difference to the higher 
efficiency of the U-spec filter, as the spectral profiles are similar. 

We then calculate the zero-point offset for the U-Bessel filter 
using the equation: 

U B es.sdzp] = Usessel ~ UBes.sc-1 ~ k V X - C.t.(U - B) 

where Usessei and usessei are the correct and observed magni- 
tudes, kjj is the extinction term for the U-Bessel filter, X the 
airmass of the observation and c.t. the colour term. To have 
consistency between the different observations made with differ- 
ent integration times, we have adopted everything to f,„, = 1 s. 
During the commissioning of the LBC-blue the extinction term 
for the U-Bessel filter was measured to be ky = -0.48 ± 0.02 
and the colour term c.t. = 0.0363 Applying the Usesseb usessei, 
and X values of the standard stars we derive: Usesseiizp] - 
(26.012 ±0.014) -2.5 log ^l. 

To calculate the U-spec zero-point offset we shift the value of 
the U-Bessel offset by the mean measured magnitude difference 
of the standard stars. This is the equivalent of assuming that the 
standard stars have the same magnitudes in the U-spec and U- 
Bessel filters. The central wavelengths and widths of the two 
filters are very similar, so such an assumption does not affect the 
result in great extent. We find: U spec [zp] = 26.073 ± 0.012. 

To calculate the zero-point offset of the final image, where 
as a result of rescaling of the individual images and combining 
them the connection to the original gains has been lost, we use 
the number mentioned in the previous paragraph to derive the 
magnitudes of Lockman Hole sources using the raw images. For 
this purpose we selected 43 non saturated sources with almost 
gaussian profiles which are observed with the second chip of the 
mosaic, as the standard stars. We use these 43 sources to mea- 
sure the zero-point magnitude of the final combined image. We 
derive zp/mai - 32.041 ± 0.013 and do not find any evidence 
for a gradient in any direction of the image. The magnitudes 

1 The LBC commissioning report can be found at 
http : //lbc . oa-roma . inaf.it/coitimissioning/ 
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Fig. 3. Final astrometry check using the USNO-A2 catalogue. The left panel shows the RA and DEC differences between the LBT 
and USNO counterparts. The (gaussian) fits to the respective histograms have cr = 0.47 arcsec and <x = 0.59 arcsec for RA and DEC 
and the means are at -0.06 arcsec and -0. 13 arcsec respectively. The right panels show the RA and DEC differences across the entire 
LBT image, where we do not observe any systematic trends. 
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Table 2. Flux calibration, quality and source extraction informa- 
tion of the images. 



of th e standard stars are given in the Landolt photometric sys- 
tem (lLandoltl [l992h which is based on Vega magnitudes. We 
use U-Bessel AB correction calculated during the commission- 
ing time (0.87), so our final zero-point offset for the U band is: 
U specAB [zp]= 32.911 ±0.013. 

Detailed information on the three (U-B-V) final images can 
be found in Tab. [2] 

5. Results 

5.1. Source catalogues 

The source detection has been done independently in each of 
the U-B-V images using sextractor. Sources are identified as re- 
gions where 12 or more adjacent pixels have values above 1.2 
times the local background rms. The algorithm first subtracts the 
background which is fitted by segmenting the image with a grid. 
If the grid is too fine, a fraction of the flux of the sources will 
be subtracted as background and this will be more severe for ex- 
tended sources. On the other hand, a very large mesh will fail 
to subtract the background near very bright objects, where stray 
light contaminates the image, so the source extraction will fail in 
these areas. To overcome these issues we ran sextractor in two 
steps: first we used a very fine grid (with a 5 x 5 pixel mesh) 



to subtract the background and created a "source detection" im- 
age. We then re-run sextractor in dual mode using this image to 
detect the sources but measure their fluxes from the original im- 
age, where the background is subtracted using a 10 x 10 pixel 
mesh. This method has the drawback that the apertures where 
we measure the flux are too tight for bright sources, which ap- 
pear more extended in the images and as a consequence we are 
losing a fraction of their flux (see Fig. [5] upper panel). Therefore 
we run sextractor once more in "single mode" (with a 10 x 10 
pixel mesh for the background; see also Fig. [5] lower panel) and 
replace the sources with magnitudes brighter than 22 of the orig- 
inal catalogue with those extracted in "single mode". 

Finally, in order to avoid spurious detections, we remove 
from our catalogues sources whose isophotal flux errors are 
larger than the fluxes and therefore do not have reliable photom- 
etry and sources whose FWHM is less than 90% of the seeing 
of each image (1.06, 0.94, and 1.03 arcsec for the U, B, and V 
images respectively) and are related with imaging artifacts. We 
also optically inspect the remaining sources and remove obvi- 
ous false detections related with bad pixels, dust on the CCDs, 
bleeding trails etc as well as saturated sources. The final U, B, 
and V catalogues contain 51500, 76071, and 68278 sources re- 
spectively. 

In order to estimate the detection limits of our catalogues we 
plot the flux error against the flux of each detected source (see 
Fig. |6j. We do this because we used a more complex selection 
algorithm to extract sources than a simple signal-to-noise cut. 
The dashed lines in Fig. |6]represent signal-to-noise ratios of 1, 
2, and 3 from left to right and the red line the (empirical) "detec- 
tion limit". We notice that the faintest sources tend to be closer 
to this limit and this is a result of them being point-like. The 
resulting detection threshold magnitudes and signal-to-noise ra- 
tios (see Tab. |2]i are not detection limits in the sense that sources 
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Fig. 4. Final U-band image and corresponding exposure map. The image covers both the XMM and the VLA fields, which are 
marked with lO'-radii circular regios, close to their typical widths. The highest exposure is in the overlapping region between them. 



with fluxes (or SNRs) above these limits are detected, they are 
indicative of the sensitivity of the survey representing the low- 
est flux (and respective SNR) of securely detected sources. They 
also provide no information on the completeness of the survey 
at a given flux (or SNR), nor an estimation of the chance of a 
spurious detection. Such an analysis is described in 35.31 

5.2. Colours 

To create colour catalogues of the various sources detected in 
the U, B, and V images, one could simply cross-correlate the 
three source catalogues described in the previous paragraphs and 
compare the fluxes in the different bands. This however would 
introduce an uncertainty on the choice of the best counterpart 
and moreover the deblending efficiency of sextractor varies be- 
tween the different images, so a source in one catalogue might 
be blended with a close pair in another or vice-versa. Therefore 
we chose to select one image to extract the sources and then 
measure their fluxes using the other images in dual mode. 

We make the source detection on a combined image of the 
three bands, the ^ 2 -image. The PSFs of the three images we 
combine do not have significant differences; the worse PSF (U- 
band) is only 6% larger than the best (B-band), therefore we do 
not lose in quality when combining the images as compared 
to using the b est PSF image and we gain in S/ N. We follow 
the recipe of Szalav, Connolly & Szokolv ( 19991) to create the 
^- 2 -image: after carefully removing any residual background of 



each image (using the "-BACKGROUND" checkimage option 
of sextractor) we fit the off-source pixel histogram with a gaus- 
sian, checking that the noise profile is indeed gaussian. We then 
scale the three images according to their noise amplitudes and 
we create the ^ 2 -image, which is the square root of the sum of 
the squares of the individual pixel values. We then extract the 
sources from the combined image using the method described 
in the previous paragraph and measure the fluxes in the indi- 
vidual U-B-V images. Again, we consider a detection real if its 
FWHM is > 90% of the PSF FWHM of the best individual im- 
age (0.85 arcsec). The final colour catalogue contains 88429 with 
detections in all three bands. 



A colour-colour (U - B vs. B - V) diagram of the Lockman 
Hole sources is shown in Fig. [7] The greyscale represents the 
density of sources detected in all three bands and the black 
lines mark the selection area of z > 3 objects. We also cal- 
culate the colours of diff e rent g alaxy SED templates from 
I Coleman. Wu & Weedmanl dl980l) and a QSO template from 
ICristiani & Viol dl990t) . The galaxy templates are extrapolated 
to the Lyman break (91 1.25 A) and are zeroed thereafter. The 
Lyman break meets the blue end of the U filter at z = 2.5 and 
this is the highest redshift where the colour tracks are reliable 
(solid lines). The dotted lines (z > 2.5) are shown as an approxi- 
mation of the colours of high redshift galaxies. 
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Fig. 6. Flux error versus flux diagrams for source detections in the U, B, and V filters. The zero-point magnitudes for the three 
filters are 32.911 + 0.013, 33.830 + 0.014, and 34.110 + 0.011 respectively. The dashed lines represent signal-to-noise ratios of 1, 
2, and 3, and the red solid lines represent the (empirical) limit of the source extracting algorithm used (4.5 <x). The horizontal and 
vertical solid lines mark the minimum flux and the respective minimum flux error of the faintest sources extracted. The dotted and 
dash-dotted lines mark the 50% efficiency limit and the limit of the number count analysis respectively. 



5.3. Number counts 

In order to derive the differential number counts of extragalactic 
sources in the U, B, and V bands, we select a region in the centre 
of the field with uniform exposure within a good approximation. 
This region has a size of 14 x 13.3 arcmin and is located in the 
area where the XMM and VLA observations overlap. 

The first step in calculating number counts is to estimate the 
source extracting efficiency at a given magnitude. The way to 
do it is to create an image with artificial sources of known mag- 
nitudes and to apply the same source extracting procedure as 
applied to the image, and measure the fraction of the sources re- 
covered. We use the artdata package in IRAF to create lists of 
artificial sources. They contain sources with a uniform spatial 
distribution and magnitudes ranging from 16 to 29 following a 
power-law distribution with a power of 0.5. The surface bright- 
ness profiles are exponential discs (resembling spirals) and r~'/ 4 
discs, resembling ellipticals. The fr action of elliptical g alaxies in 
the random catalogues is 20% (see Ivan den Berghil200l . Here, 
we caution that adding a large number of artificial sources in the 
image might change its crowding properties, however we need 
a large sample of sources for reliable statistics. To avoid con- 
fusion, we create a list of 100000 sources and split it to 100 
1000-source samples. 

We plant these sources into the cutouts of the final images 
and apply the same source extracting algorithm we used to cre- 
ate the source catalogues. We then measure the fraction of the 
artificial catalogue we retrieve, hence the efficiency of the source 
detecting method at any given magnitude and average the results 
of the 100 subsamples. Increasing the number of sources of each 
subsample we get similar results up to the point where the num- 
ber of sources is comparable to the number of "real" sources in 
the region (~ 10000). The results for all three bands are shown in 
Fig. [8] This method has the drawback that the artificial sources 
are mixed with real sources, and so there is no way of knowing 
whether a detected source is real or an artifact. The surface den- 
sity of sources with magnitude (at any band) < 27 is close to 
2 x 10 5 deg~ 2 , which means that there is a ~ 10% probability that 
a real source is within 1.5 arcsec of a random position. 

To measure the spurious source detection rate we measure 
the off-source noise of the science images and check that the 
noise profile is gaussian. We then create gaussian noise maps 
of the same amplitude and insert the artificial sources there. In 
this case it is desirable to reproduce the crowding of the original 
field, so we include a large number of artificial sources (30000), 



which is the number of sources with magnitude < 29 we ex- 
pect in this field. The source extraction output to these composite 
images provides the information of the spurious detection rate, 
plotted with the red lines in Fig. [8] We can see that the num- 
ber of spurious sources is negligible below 26.5 mag and starts 
becoming important above 28.0 mag, where the efficiency drops 
to practically unusable values. From these diagrams we can also 
derive the magnitude where the detecting efficiency drops below 
0.5, which is a meaningful measure of the detection threshold of 
the image. This threshold is 26.7mag(AB), 26.3mag(AB), and 
26.3 mag(AB) for the U, B, and V bands respectively. 

Figures l9lfTT1 show the surface density distributions for the 
three bands observed. Data points of other studies found in the 
literature are also plotted. We have binned the magnitudes of 
the observed sources in bins of 0.5 mag. For each bin we cor- 
rected the source counts using the efficiency and spurious de- 
tection information. As we are interested in galaxy counts, we 
need a selection mechanism for stellar sources, and as such we 
use the "stellarity index" of sextractor. This estimate works well 
for bright sources, but because fainter galaxies can appear point- 
like it fails for larger magnitudes. As a limiting magnitude we 
choose 21.5(AB). Below this limit the stellar counts are any- 
way negligible with respect to th e number of galaxies (see e.g. 
iJarrett. Dickman & Herbstl [1994). The error bars take into ac- 
count Poisson uncertainties of the uncorrected counts, efficiency 
uncertainties and cosmic variance. For the latter we use the com- 
putational tool of iTrenti & SfiaveiiH d2008l) . which compares the 
two-point correlation function of dark matter with the volume of 
the survey. As a typical redshift for our survey we use z — 1 (see 
also {36. Il l, though using a different value in the range 0.5 - 2.5 
does not change the result in great extent. We find that the cos- 
mic variance uncertainty is important for bright magnitudes (typ- 
ically < 25.5) where the number of intrinsic objects is relatively 
small. For fainter magnitudes the efficiency uncertainties dom- 
inate. As an estimation of those we choose the mean efficiency 
difference between the bin in question and its neighboring bins. 
This way we account for the effects of binning, in other words 
the different efficiencies the magnitudes within each bin have. 

The surface density data can be seen in Tab. [4] We com- 
pare them with the results of various other studies (Figures [9] 
[TUl [TT1 Tab. [3j and we are in good agreement. The small di- 
agrams of Figures [9j [10] andQT]plot the depth reached by each 
survey presented with respect to its covered area (the various sur- 
veys used to create these figures are presented in Table [3j. The 
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Fig. 8. Detection efficiency histograms of the source detection algorithm used, derived from simulations. The blue lines represent 
the smoothed efficiencies used to correct the surface density distributions. The red histograms are the estimates of the fraction of 
detections which are spurious, using the artificial noise images. 



LBT survey presented here is the deepest one in the U band ever 
conducted in such a large area and among the deepest in the B 
and V bands. Significantly lower limits have been achieved only 
with the Hubble Space Telescope in pencil-beam surveys (HDF- 
N and H DF-S), and these are high ly sensitive to cosmic vari- 
ance (see Sommerv ille et all 120041) . Comparing or results with 
those from surveys of similar widths (made with the LBT and 
the Subaru telescope) we find very good agreement. 



6. Discussion 



6.1. Number counts 



The models com piled by Metcalfe et al.l ([1996) and 
iMetcalfe et al.l d200ll) (normalized to 18 mag using all data 
available in the bibliography) are plotted against our mea- 
surements for the U and B bands in Fig. FPU We plot here 
three of the models presented in Metc alfe et al.l (1 19961) and 
IMetcalfe et al.l (12001 1) . The short-dashed lines represent the 
pure luminosity evolution model, the long-dashed line the same 
model with the inclusion of a population of star-forming dwarf 
galaxies, which is the best-fit model in Metc alfe et al.l 01996) 
and the solid line is the same pure luminosity evolution model 
with a modification of the faint-end slope of the luminosity 
functions of late-type spirals (a = — 1 .75 instead of a = - 1 .5), 
used to fit the multi-colour data of IMetcalfe et all (120011) . We 
find very good agreement with the a = -1.75 model in the 
B-band, although the U-band counts are under-predicted by all 
models. However, the faint-end slope of the U-band counts does 
seem to support a steepening of the faint-end slope of the LF. 
Barr o et all (120091) have shown that the slope of the number 
count distribution assymptotically reaches -0A(a +1), where a 
is the faint-end slope of the luminosty function if parametrized 
by a Schechter function. Measuring the slopes of the number 
counts using the five faintest points of each distribution, we 
calculate the faint-end slopes of the respective luminosity 
functions: a v = -1.733 + 0.018, a B = -1.748 + 0.006, 
ay - -1.507 + 0.018. We note that the assumed steep LF 
faint end slope is in good agreement with the number count 
distributions of the U and B bands, whereas the V band points 
to a LF with a = -1.5 (see IMetcalfe et al.Ll200lb . 

At this point it is useful to have a notion about the type of 
galaxies that are best represented in our sample and their red- 
shifts. A valuable tool in this direction is the colour distribution; 
Fig. [7] shows the colour plot of the Lockman Hole sources with 
tracks of templates of different types of galaxies. We note that 



the spiral (Sbc-Scd) and irregular tracks lie closer to the bulk 
of observed colours. The metallicity and extinction properties 
of a galaxy can have a severe effect in its optical colours. For 
that reason we reproduce Fig. [7] with a set of SED templates 
which have varying stellar ages, met allicities and extinctions , 
calculated with the GISSEL98 code dBruzual & Charlotill993"h . 
The results are presented in Fig.[12j where the colour tracks are 
colour-coded with respect to the redshift. The distribution of 
sources is well reproduced and we can see that the redshift range 
most represented is z — 1 .5 + 1 . Moreover, the bulk of the colour 
distribution is represented by spiral and irregular tracks, while 
the ellipticals account for colours redder than B —V — 1 . 

The galaxy types and redshift probed by our s urvey are com- 
patible with the "steep faint end slope" model of IMetcalfe et al.l 
(2001) and the slopes are also in good agreement. Therefore, 
there is no need to invoke a dwarf galaxy population to assist 
the sources which cause the steepening of the faint end slope, 
in order to reproduce the B-band data. The U-band counts on 
the other hand are underestimated by the "steep faint end slope" 
model, although the slope itself agrees. In this case a dwarf 
galaxy population would assist in incrasing the U-band number 
counts. It would be however challenging, as such a population 
is required to affect the U-band leaving the B-band unchanged. 
The Ly-ff line falls into the U wavelength range at a redshift 
of z — 2, so a population if highly ionized Ly-a emitters is a 
good candidate. However, at z = 2.5 the blue filter would also be 
affected, which leaves a narrow redshift window for this hypo- 
thetical population. An implication of this scenario is a sizeable 
decrease in t he star formation rate between z — 2 and z = 2.5. 
iReddv et al. (2008) find an increase in the star formation rate 
density between z ~ 3 and z ~ 2, which is reflected in the UV 
luminosity density, and more specifically in the number density 
of faint (M AB (1700 A) > -21) UV-emitting galaxies. 

A steepening of the faint end slope of the (B-band) luminos- 
ity function is already ev ident since the fir st computation of its 
values at redshifts z > dLillv et al.L [T995). These authors find 
that the slope increases with redshift (out to z = 1.3) and that 
this is an effect caused by galaxies with blue optical colours, 
while the LF of red galaxies shows minimal change in its fit- 
ted parameters. This result is back e d up by more recent s tudies 
" . .2004 lllbert et all 120051: lArnouts et all 120051: 
2006:" 1Prescott. Baldrv & Jamesl l2009TT and the 



(Gabasch et al 
Willm er et al 



general trend is that the not only blue galaxies' LFs evolve more 
with redshift than those of redder colours, their faint end slopes 
are steeper as well. In cases where the LFs are computed with re- 
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Table 3. Number counts data found in the literature. The symbols noted are used in Figuresl9lfTT1 A "d- 
symbolizes digitization of photographic plates. 



prefix before the instrument 



spect to the galaxy type dllbert et all 120061; IZucca et all [2006), 
little (if any) evolution of the faint end slope is found for each 
galaxy type, while th e slope is different f or each type, the sttepest 
beeing in ir r egular s dZucca et all 120061) or blue-bulge galaxies 
dllbert et all 120061) . There is however significant change in the 
normalization and the value of M* (the characteristic Schechter 
luminosity), which is interpreted as an increase in the fraction of 
irregular and late-type galaxies with redshift. The steep faint-end 
slope we find in the U and B bands ( a -1.75) agrees w ith the 
values fitted for irregular galaxies bv lZucca et alj d2006l) and is 
even a bit too flat compared to the value assumed bv lllbert et alj 
(2005:?) for the blue-bulge population (their a = -2.0). Given 
that in this survey the dominant population, especially at faint 



magnitudes, is spirals and irregulars at non-local redshifts, we 
support these steep faint-end slopes. The sources responsible for 
the steep slopes are activly star forming and are goo d candidates 
for the "blue dwarf" population. lDriver et alj d 19951) assume that 
this population consists of sourc es with late-type and irregular 
morphologies; Ilbe rt et aT] d2006l) state that the "blue bulge" pop- 
ulation could be a population of actively star-forming galaxies, 
where the starburst region has bulge-like morphologi cal char- 
acteris tics, like the "blue spheroid" galaxy sample of llm et all 
(120011) . 

An issue that still needs to be addressed is the flattening of 
the number counts slope in the V-band. A mechanism that affects 
the faint-end slope is supernova feedback, which is caused by 
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Table 4. Differential number counts (in 0.5 mag bins), uncertainties and source detections efficiencies for each magnitude bin of 
galaxies in the U, B, and V bands 
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Fig. 13. Measured su rface densities in the U and B band (red points), plotted with evolution models from Metc alfe et al.l d 19961) and 
Metc alfe et al.l d200ll) . The short-dashed lines represent the pure luminosity evolution model, the long-dashed lines the same model 
with the inclusion of a population of star-forming dwarf galaxies, and the solid lines the same pure luminosity evolution model with 
a steep faint-end slope of the luminosity function. The latter is in very good agreement with the B-band counts, while the U-band 
counts are under-predicted by all models. 



the he ating of the interstellar medium through supernova expo- 
sions. Nagashima et al. (2005) have modelled galaxy formation 
taking this effect into account. Their predictions for the number 
counts using strong or weak feedback (parametrized by the time- 
scale in which supernova explosions reheat the cold interstellar 
gas) differ in the faint slo pe with minimal impact on the normal- 
ization (see Figure 18 in Nagashimaetall 12005b. Fig.[T4lplots 
the B and V-band number counts predictions of Nagas hima et all 



(2005) with our data-points; the solid and dashed lines refer to 
strong and weak feedback respectively. While both predictions 
seem to overestimate the observed number counts at faint fluxes, 
the faint end slope of the weak feedback prediction is in good 
agreement with the data for the B-band, while the V-band slope 
is better interpreted with the strong feedback model. A possi- 
ble explanation is that the V-band probes the same rest frame 
wavelength at higher redshift. If we assume that star formation 
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Fig. 5. Apertures of extracted sources in the U image. The up- 
per panel shows the "dual mode" results and the lower the "sin- 
gle mode" (see text) Note that the detection efficiency near the 
bright object as, the deblending efficiency and the apertures of 
faint sources are much better in "dual mode", but the apertures 
of bright objects are more reasonable in "single mode". 



is the dominant mechanism producing near-infrared light (where 
the rest-frame B and V bands are at redshift z > 1.5) the V- 
band probes higher redshifts than the B-band. There is evidence 
that the UV luminosity func tion has a steeper f aint end slope 
at z ~ 2 (a = -1.8 8 + 0.27.lReddv et all 120081) than at z ~ 3 
(a = -1.60 + 0.13. ISteidel et all 1 1 999jT ln this case, starburst 
feedback would be stronger at higher redshift, in line with our 
data. So, enhanced SFR at z > 1.5 could cause the flattening of 
the V-band faint slope. 

6.2. Colour selection 

To be able to test evolutionary models of galaxies in a greater 
extent one needs to have information of the redshifts of the var- 
ious objects found in a "blind" survey. However, even with the 
largest telescopes available it is practically impossible to have 
complete samples beyond R ~ 24 and use the full capacity of 
photometric surveys. Moreover, the selection of targets for spec- 
troscopy at such faint limits is hard because their redshift range 
is so large that it makes it impossible to get meaningful spec- 
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Fig. 7. Optical colours of the Lockman hole sources in the AB 
system. THe greyscale represents the density of sources detected 
in all three bands are plotted. The black line marks the selection 
area of U-dropout sources with redshift z > 3 (see Fig. [12] and 
§ 16. 2b . Dotted are the tracks above z - 2.5, where the SEDs are 
not well sampled. The coloured lines represent tracks of different 
SED templates: blue line for ellipticals, cyan for irregular, green 
for Sbc, yellow for Scd, and red for QSO. 
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Fig. 9. U-band surface density distribution corrected for effi- 
ciency and spurious detections. The data points of this study 
(seen in Tab. @|i are plotted in red dots, while various symbols 
represent the results of other surveys found in the literature (see 
Tab. [3]), corrected to AB magnitudes. The inner plot shows the 
depth of each survey (in terms of secure number count measure- 
ments) with respect to its width. 
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Fig. 14. Measured surfa ce densities in the B and V band (red points), plotted with number counts from the simulations of 
Nagashi ma et al.l ((2005). The solid and dashed lines represent the strong and weak supernova feedback cases respectively. The 
faint end slope of the weak feedback model is in better agreement with the B-band data, while the the strong feedback model is in 
better agreement with the V-band data. 
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tra without pre-selecting the targets according to their redshift 
range. A way to overcome this barrier is to use the photometric 
redshift technique, where the SED of each source is compared 
with known SED templates to derive an estimate of the redshift. 
Although the accuracy of this method is limited so it cannot be 
used for e.g. spatial clustering studies it can be very useful in 
deriving luminosities or selecting objects in different redshift 
ranges. Given the extensive spectral coverage of the Lockman 
Hole it is possible to calculate photometric redshifts for a large 



number of galaxies. Details about the Lockman Hole photo-z 
survey will be given in a subsequent paper. 

The drawback of the photometric redshift technique is that it 
requires the detection of the source in a large number of bands 
spanning from the near ultraviolet to the infrared. It is how- 
ever possible to select so urces within a redsh ift range using the 
"dropout" technique (e.g. lSteidel et aUl2003l) . This technique is 
used to detect the Lyman break in the spectra of galaxies when 
it is redshifted between two of the observed bands. In practice 
the colour-colour diagram is used to select the sources. Based on 
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Fig. 12. The same plot as Fig. [7] but with with SE Ds created 
with the GISSEL98 code dBruzual & CharlotLll993"l) with differ- 
ent m etallicity and extinction properties (see iBolzonella et al.L 
2000). We use elliptical, irregular and spiral tracks, colour-coded 
with redshift ranges. The bulk of the distribution is reproduced 
with spiral and irregular tracks at redshift 0.5 < z < 2.5. The 
z > 3 region is also marked. 



the colour tracks of Fig. [12] we set the selection limits of sources 
with z > 3 to: 

U-B>0.9 

U -B> B-V + 1.2 

B-V <13 

Using these limits, we find 2152 sources with redshift z > 
3 in the whole Lockman Hole regi on (925 arcm i n 2 ), w hich 
gives a number density of 8375 deg" 2 . ISteidel et al.1 d2003l) find 
6176 deg~ 2 U-dropout z ~ 3 sources using different selection 
bands (U„-G vs. G-'R). Using the near-infrared coverage of the 
Lockman Hole we could also s e lect s ources according to their 
B-z-K colours (see iDaddi et all 12004), having 1.4 < z < 2.5, 
or according to their i-z colour, having 1.4 < z < 2.5 (see 
IVanzella et all 120051) . Being able to select sources in distinct 
redshift ranges provides a valuable tool to further test models 
of galaxy evolution. 

Selecting sources by their colours can also provide sam- 
ples of different kinds of objects. Compton thick AGN are ac- 
tive galaxies with dense environments (TVe > 10 24 cirT 2 ) so 
that they block even hard X-ray radiation and are not detected 
even in the deepest X-ray surveys. They would provide valu- 
able information in evolution studies, as they represent a dis- 
tinctive phase of a galaxies lifetime and they are the "miss- 
ing lin k" in population synthesis models of the X-ray back- 
ground dGilli. C omastri & HasingerL l2007l) . The most promising 
methods of detecting Compton thick AGN involves comparing 
the optical and inf r ared fluxes o f sour c es dDonlev et all 12 007 ; 
IDaddi et al.L 120071; iFiore et al.L I2008L iGeorgantop oulos et all 



2008). The multi wavelength coverage of the Lockman Hole in 
combination with the deep X-ray observations are ideal for this 
kind of study. 

6.3. Follow-up 

One important contribution of this study is that it provides a large 
number of newly detected extragalactic objects to be further ob- 
served in follow-up campaigns. A number of sources has already 
been spectroscopically identified, and the y have been s e lected 
from the X-ray campaigns with ROSAT ( Sch midt et ail 1 1998L 
Lehmann et ail I2000L 1200 ll) and XMM-Newton dMateos et all 



2005). With multi-object spectrographs we are now able to con- 



duct spectroscopy to a large number of optical sources. The 
LBT is already equ ipped with a near-in frared multi-slit spectro- 
graph (LUCIFER; iMandel et ail 120071) which will start opera- 
tion within 2009 an d the optical multi-slit spectrograph (MODS; 
Pog ge et ail 120061) is expected to be operational in 2010. The 
key scientific goals of these instruments is to conduct spec- 
troscopy at cosmologically interesting redshifts. To be able to 
select targets for these instruments we need a deep optical sur- 
vey and a colour selection scheme similar to what described in 
the previous section. 



7. Summary and conclusions 

In this paper we present the deep imaging campaign of the 
Lockman Hole using the LBT. The Lockman Hole is an excellent 
region for deep multi-wavelength observations given the mini- 
mal galactic absorption. Here we report details of the U, B, and 
V-band observation and the data reduction strategy. Our imag- 
ing area covers 925 arcmin 2 in a very well sampled region of 
the Lockman Hole, with deep X-ray, infrared, and radio cover- 
age. We have reached depths of 26.7, 26.3, and 26.3 mag(AB) 
in the U, B, and V band respectively, in terms of 50% source de- 
tection efficiency, and have extracted a large number of sources 
(~ 85000) an all three bands. 

The number counts distributions are used to test galaxy evo- 
lution models and and simulations. We find evidence of steep- 
ening of the faint-end slope of the luminosity function in the U 
and B bands, which can explain the B number count without the 
need of a dwarf galaxy population. However the U counts are 
under-predicted with this model and an enhancement of the star 
formation rate at z = 1.5 -2.5 is needed to explain them. A flatter 
faint end slope observed in the V-band case could be the result 
of supernova feedback. 

This survey is part of an effort to conduct deep observations 
of the Lockman Hole in different bands ranging from the infrared 
to the X-rays. This will help us select different source classes for 
further study and in addition to planned spectroscopic observa- 
tions create a large database for extragalactic studies. 
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